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Abstract 

The production of J/ip at large transverse momenta (p± > Mju,) in pp collisions is 
considered by including the mechanism of fragmentation. Both contributions of fragmen- 
tation to J/0 and of fragmentation to x states followed by radiative decay to J/ip are 
taken into account. The latter is found to be dominant and larger than direct production. 
The overall theoretical estimate is shown to be nearly consistent with the experimental 
observation. 
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Figure 1: One of the diagrams for 
the gluon fragmentation function at 
the scale \x = 2m c . 



Figure 2: One of the diagrams for 
the charm fragmentation function, 
at the scale \i = 3m c . 



The study of the properties of the bound states of heavy quarks plays a central role in 
the understanding of Quantum Chromodynamics in that it stands on the very border between 
the perturbative and non-perturbative domain. In particular it is of key importance to have 
accurate estimates of the production cross sections at large tranverse momenta for precision 
tests of the theory and possible evidence of new phenomena. 

So far there has been an intensive experimental study of the cc IS vector bound state, 
namely J/0, in hadron collisions both at UAl and CDF 0. The results have been compared 
with theoretical calculations which take into account two different mechanisms for J/ijj 
production: direct charmonium production, including the contribution from the \ states, i.e. 

gg^jUg "'^ ^ X9 19 ^ XI ^ 

and the production resulting from B mesons decay 

pp — > bX 

U B -»■ J/ip X 



(2) 



A more recent version of this calculation, which makes use of the NLO prediction Q for b 
production, is presented in ref. p. 

These calculations are in disagreement with the results form CDF, the J/ip rate observed 
being actually higher, by a factor of two or more, than the predicted one 0, |5], |6|]. 

It has however recently been pointed out by E. Braaten and T.C. Yuan p| that at large p± 
an additional production mechanism comes into play, namely the fragmentation of a gluon or 
a charm quark into a charmonium state. While being of higher order with respect to direct 
production by a power of the running coupling constant a s , this mechanism becomes dominant 
at large p± because of a factor 0{p 2 L /m 2 ) which overcomes the extra power of a s . The frag- 
mentation functions describing these processes can be calculated perturbatively. Indeed it has 
been argued in J7| and subsequently shown at LO in |J that J/ip production via fragmentation 
will overcome the direct one (i.e. gg — > J/ip g) at p± ~ 6-8 GeV. A similar exercise for the 
X production, when the total fragmentation probability / (see below) times the gg — ► gg 
cross sections is compared to the direct production gg —>■ xg, reveals that fragmentation should 
dominate for p± already at ~ 2 GeV. Since this result is at the limit of validity of the fragmen- 
tation function approach, we can however still expect that the fragmentation mechanism will 
dominate over the direct one at p± values as low as 5-6 GeV. 

In this Letter we apply these ideas to a quantitative determination of the J/ip production 
rate in hadron collisions, taking also into account the production via fragmentation processes 
of the x states and subsequent radiative decays to J/ip. 
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Figure 3: One of the diagrams for 
the gluon fragmentation function to 
the x states, at the scale \i = 2m c . 



Figure 4: One of the "perturba- 
tive" contributions to the induced 
gluon fragmentation function, at the 
scale (j, = 4m c . 



To this aim the following fragmentation functions play a major role: the gluon fragmentation 
to J/ip J7j, DgW (see fig. [I]); the charm (or anticharm) fragmentation to J/ij) ||, D J J^ (see 
fig. ||); the charm fragmentation to x states |1J, D*; and finally the gluon fragmentation to x 
states 0, D* (see fig. |). 



They have been all calculated by perturbative techniques at an initial scale of the order 
of the mass of the J/ip. Of course in the evaluation of the actual cross sections they must be 
evolved to the appropriate scaleQ, and one gets to the usual expression 

da(pp - J/VKpx) +X)=J2[ 1 *z MPP - Hpjz) + X, fi)Df(z, /x) (3) 

for the J ftp production, the sum running over g, c and c. A similar formula does hold for x 
production?]. The cross section on the right hand side corresponds to the inclusive production 
of the parton i, convoluted with the appropriate structure functions (throughout this work we'll 
use the MRS-DO set), and summed over all relevant parton-parton scattering processes. \i is 
the factorization scale, which we will take of order /x = \jp\ + Mj^. 

The evolution of the fragmentation functions given above obeys the usual Altarelli-Parisi 
(AP) equations 

H^-Df{z^) = E /' - P^/Viti Df{y,n) (4) 



Furthermore it has been pointed out in ref. (T^j that when one considers the whole set of the AP 
equations, with the appropriate mixings taken into account, the evolution of the D 3 J^ will induce 
a gluon fragmentation function through the splitting g cc and subsequent fragmentation of 
one of the quarks into a Jfip (see fig. |j). This process is of order but, being enhanced by a 
factor log(/x/Mj/^), will dominate over the the contribution from D J J^ at large p±. 

We present our results first using the leading order (LO) expressions for the partonic cross 



sections and then, to reduce the theoretical uncertainty, by taking into account the full NLO JTjJ ] 
information on the partonic scattering processes. 

We plot in fig. | the LO cross sections, differential in p± and integrated over the rj < |0.5| 
range, for producing a J/0 via fragmentation, either directly or after radiative decay of a x state. 
The values of the various parameters entering into the calculation are listed in Table p], and 

1 We thank Paolo Nason for having provided us with the Fortran code for the numerical evolution of the 
fragmentation functions. 

2 Note that, when performing the convolution with the xo,2 fragmentation functions, due to the singularity 
that these functions display at x — 1 an appropriate subtraction precedure, already introduced in ref. []l2|| , must 
be used to ensure the convergence of the numerical integration. 
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a s (2m c ) 


0.26 






1.5 GeV 


Ref. [0 


\Ro(o)\ 2 


(0.8 GeV) 3 




Hi 


15 MeV 




H' s {m c ) 


o ivre v 


Ref. [11 




BR(xo 




0.007 




BR(xi -> 




0.027 


Ref. [15 




BR(X2 - 




0.014 




BR(J/ty> -» 


0.0597 




Initial scale for 

charm fragm. evolution (fJ% n i) 


3m c 


Ref. § 


Initial scale for 

gluon fragm. evolution (/if m -) 


2m c 


Ref. § 


Initial scale for 

induced gluon evol. (A4ni ) 


4m c 


Ref. [16 





Table 1: Summary of the parameters. We have used for 
the fragmentation functions the same values given in the 
cited references, and the branching ratios as quoted by the 
Particle Data Group. 



we are using fi = fi for the factorization/renormalization (f/r) scales. The curves labeled by x 
are due to gluon fragmentation only. We have not included c — > \ fragmentation contributions 
since, from the total fragmentation probabilities listed in |ID[ , they can be predicted to lie two 
orders of magnitude below the c — > J/ip curve and be therefore surely negligible. 

From inspection of fig. |5] the contributions from xi an d Xi states can be clearly seen to 
dominate all over the p± range considered. 

Next we compare, in fig. |6], the results obtained for the dominant Xi + X2 contribution in the 
LO approach with those obtained by inserting also the next-to-leading (NLO) partonic cross 
sections, to order a 3 , with a s evaluated to two loop accuracy. We remark that the latter does 
not result in a full NLO calculation as the initial state evaluation and the evolution kernels 
of the fragmentation functions are still at leading order. We have explicitly checked that the 
inclusion of NLO evolution kernels does not change appreciably our results. Figure ^| clearly 
shows that the higher order terms enhance the cross section by a factor about 1.5. This is 
consistent with previous studies of higher order corrections in heavy quark 0, W% and inclusive 



jets fi~4] , |l~8f production in hadron collisions. The effect of variations of the f/r scales fi between 
0.6yU and 2/z is also shown. As expected, the inclusion of the NLO terms reduces the sensitivity 
to scale variations. 

One more theoretical uncertainty has to be considered, namely the freedom to choose the 
initial scale for the fragmentation functions evolution around Mj^. Figure shows the effect 
of varying fif ni in the range 2-4 GeV. 

Finally we show, in figure 181 and in table 0, our final prediction for J/ip production by 
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direct (l9| + B decays f§ 


fragmentation 


total 


p± 


lower 


upper 


lower 


upper 


lower 


upper 


5.5 


.52 x 10° 


.10 x 10 1 


.12 x 10 1 


.28 x 10 1 


.17 x 10 1 


.39 x 10 1 


7.5 


.17 x 10° 


.36 x 10° 


.30 x 10° 


.75 x 10° 


.46 x 10° 


.11 x 10 1 


9.5 


.68 x KT 1 


.14 x 10° 


.98 x 10" 1 


.25 x 10° 


.17 x 10° 


.39 x 10° 


11.5 


.29 x KT 1 


.67 x 10" 1 


.39 x 10" 1 


.10 x 10° 


.67 x 10" 1 


.17 x 10° 


13.5 


.14 x KT 1 


.31 x 10" 1 


.17 x lO" 1 


.47 x 10" 1 


.32 x 10" 1 


.79 x 10" 1 


15.5 


.77 x 1(T 2 


.18 x 10" 1 


.86 x 10~ 2 


.24 x 10" 1 


.16 x 10" 1 


.41 x 10" 1 


17.5 


.43 x 10" 2 


.99 x 1(T 2 


.46 x 10~ 2 


.13 x 10" 1 


.89 x 10" 2 


.23 x 10" 1 


19.5 


.25 x 1(T 2 


.62 x 10~ 2 


.26 x 10" 2 


.73 x 10~ 2 


.52 x 10~ 2 


.13 x 10" 1 



Table 2: Numerical values of the cross sections (nb/GeV) 
plotted in fig. | (BR(J/0 -> n + n~) included). 



adding the mechanism of fragmentation to the direct one 19 and to the production from B 
decays as taken from ref . || , together with the reported theoretical uncertainty. The bands are 
made by choosing the highest and lowest curve which could be obtained by varying some of the 
parameters: the factorization/renormalization scale and the value of A in the work of ref. |5|, 19 . 
again the same scale and the initial scale \i\ ni in our result. The values of the parameters used 
for each curve are given in the caption. The total result is obtained by adding togheter the two 
highest and the two lowest curves respectively. The size of the the fragmentation contribution is 
seen to be comparable with the previous estimate for the sum of the two mechanisms considered 
up to now, leading therefore to a sizeable enhancement of the predicted overall production rate, 
which we also show in the figure. 

When we finally compare with CDF data points [[J we see that they are now more compatible 
with the theoretical band. This improves sensibly the previous situation, where only by making 
very extreme choices of the parameters one could get close to the experimental findings. 

To conclude, we have considered the inclusive production of J/ip in hadron collisions in 
the framework of fragmentation functions. We have shown explicitly that the production and 
successive radiative decays of the \ states plays a dominant role. The overall theoretical 
estimate, including the contribution from B decays, is nearly consistent with the experimental 
observations. 

A more detailed analysis, and a comparison with the most recent Tevatron data, will be 
presented elsewhere pp[ . 

One of us (M.C.) would like to thank Stefano Moretti, Marzia Nardi and Roberto Tateo for 
useful suggestions. 
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Figure 5: Leading order differential cross sections due to 
various fragmentation processes: c, g, induced gluon frag- 
mentation to J/ip and gluon fragmentation to x followed by 
radiative decay to J/ip are shown. Parameters as in table H. 



Figure 6: Leading vs. next-to-leading cross section for 
producing a J/tp via fragmentation. Only the dominant xi 
and X2 contributions are included. 
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Figure 7: Theoretical uncertainty due to changes in the 
initial scale \j? ini for the fragmentation functions evolution. 



Figure 8: Theoretical prediction for J/tp production at 
Tevatron, compared to experimental data from CDF [Q]. 
Both the old result (dotted line) from ref. [||, (upper 
curve: [i = ttit/^, A = 275 MeV; lower curve n = uit, A = 
215 MeV) and the new fragmentation (dashed line) contri- 
bution (upper curve: \x = 0.6/io ; M?nj = 4 GeV; lower curve: 
H = 2^0) y^ni = 2 GeV) are included. 
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